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Abstract 

Recently, two types of simulations (forward and backward simulations) and four types of bisimulations (forward, back- 
ward, forward-backward, and backward-forward bisimulations) between fuzzy automata have been introduced. If there 
is at least one simulation/bisimulation of some of these types between the given fuzzy automata, it has been proved that 
there is the greatest simulation/bisimulation of this kind. In the present paper, for any of the above-mentioned types of 
simulations/bisimulations we provide an effective algorithm for deciding whether there is a simulation/bisimulation of 
this type between the given fuzzy automata, and for computing the greatest one, whenever it exists. The algorithms are 
based on the method developed in [J. Ignjatovic, M. Ciric, S. Bogdanovic, On the greatest solutions to certain systems 
of fuzzy relation inequalities and equations. Fuzzy Sets and Systems 161 (2010) 3081-3113], which comes down to the 
computing of the greatest post-fixed point, contained in a given fuzzy relation, of an isotone function on the lattice of 
fuzzy relations. 
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1. Introduction 

Bisimulations have been introduced by Milner |41] and Park [46] in computer science, where they have 
been used to model equivalence between various systems, as well as to reduce the number of states of these 
systems. Roughly at the same time they have been also discovered in some areas of mathematics, e.g., 
in modal logic and set theory. They are employed today in a many areas of computer science, such as 
functional languages, object-oriented languages, types, data types, domains, databases, compiler optimi- 
zations, prog ram anal ysis , verification tools, etc. For more information about bisimulations we refer to 

Jl [il IHIIIm SI BMlii . 

The most common structures on which bisimulations have been studied are labelled transition systems, 
but they have also been investigated in the context of deterministic, nondeterministic, weighted, probabil- 
istic, timed and hybrid automata. Recently, bisimulations have been discussed in the context of fuzzy auto- 
mata in Illl.ll5l4l7ll32l.l5llj64lj65ll . One can distinguish two general approaches to the concept of bisimulation 
for fuzzy automata. The first approach, which we encounter in [11, 51, 65], uses ordinary crisp relations and 
functions. Another approach, proposed in flB - ItL 33,164 ], is based on the use of fuzzy relations, which have 



been shown toprovide better results both in the state reduction and the modeling of equivalence of fuzzy 
automata. In lll6l,[l7ll the state reduction of fuzzy automata was carried out by means of certain fuzzy equi- 
valences, which are exactly bisimulation fuzzy equivalences, and in [64] it was proved that even better 
results can be achieved by using suitable fuzzy quasi-orders, which are nothing but simulation fuzzy quasi- 



orders. Moreover, it turned out in 116lll7Ll64ll that the state reduction problem for fuzzy automata is closely 
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related to the problem of finding solutions to certain systems of fuzzy relation equations. This enabled not 
only to study fuzzy automata using very powerful tools of the theory of fuzzy sets, but also, it gave a great 
contribution to the theory of fuzzy relational equations and has led to the development of the general theory 
of weakly linear fuzzy relation equations and inequalities in IsbI] . 

The same approach has been used in [15], in the study of simulations and bisimulations between fuzzy 
automata, where simulations and bisimulations have also been defined as fuzzy relations. There have been 
introduced two types of simulations (forward and backward simulations) and four t5rpes of bisimulations 
(forward, backward, forward-backward, and backward-forward bisimulations). There has been proved that 
if there is at least one simulation/bisimulation of some of these types between the given fuzzy automata, then 
there is the greatest simulation/bisimulation of this kind. However, there has not been given any effective 
algorithm for deciding whether there is a simulation/bisimulation of some of these types between the given 
fuzzy automata, and for computing the greatest one, if it exists. In [W} a theorem has been proved which 
can be used just for checking the existence of a uniform forward bisimulation (i.e., a complete and surjective 
forward bisimulation) between the given fuzzy automata. According to this theorem, there is a uniform 
forward bisimulation between fuzzy automata and S if and only if there is a special isomorphism 
between the factor fuzzy automata of ^ and S with respect to their greatest forward bisimulation fuzzy 
equivalences. Using the algorithm provided in [1^, in numerous cases we can effectively compute the 
greatest forward bisimulation fuzzy equivalences E on J?[ and F on S. Then we can construct factor 
fuzzy automata ^/E and S/F, and check whether there is an isomorphism between them that satisfies an 
additional condition stated in fl5]. But, even when we are able to effectively compute the greatest forward 
bisimulations E and F and construct the factor fuzzy automata ^/E and S/F, it maybe difficult to determine 
whether there is an isomorphism between ^/E and S/F that satisfies this additional condition. This problem 
comes down to the well-known graph isomorphism problem, which is one of the few important algorithmic 
problems whose rough computational complexity is still not known, and it is generally accepted that it 
lies between P and NP-complete if Pt^^NP (cf. [63]). Fortimately, although no worst-case pol5momial-time 
algorithm is known, testing graph isomorphism is usually not very hard in practice. 

In this paper, for any of the above mentioned types of simulations/bisimulations we provide an effective 
algorithm for deciding whether there is a simulation/bisimulation of this type between the given fuzzy auto- 
mata, and for computing the greatest one, whenever it exists. The algorithms are based on the method devel- 
oped in [,33], which comes down to the computing of the greatest post-fixed point, contained in a given fuzzy 
relation, of an isotone function on the lattice of fuzzy relations. Namely, for each type of simulations and 
bisimulations we determine the corresponding isotone and image-localized function (p on the lattice of fuzzy 
relations, as well as the corresponding initial relation ip, and the computing of the greatest simulation/bisim- 
ulation of this type we reduce to the the computing of the greatest post-fixed point of cp contained in \p. This 
is an iterative procedure by which we successively build a decreasing sequence of relations, starting from 
the relation ^ and using the function (p. If this sequence is finite, then it stabilizes and its smallest member is 
exactly the fuzzy relation which we are searching for, the greatest post-fixed point of (p contained in ip. We 
determine sufficient conditions under which this sequence is finite, when our algorithm terminates in the 
finite number of steps (cf. Theorem 15.31 1, as well as sufficient conditions under which the infimum of this 
sequence is exactly the fuzzy relation which we are searching for (cf . Theorem l5.4b . Modifying the algorithms 
for computing the greatest simulations and bisimulations we provide algorithms for computing the greatest 
crisp simulations and bisimulations between fuzzy automata (cf . Proposition l5.5l l. These algorithms always 
terminate in a finite number of steps, independently of the properties of the underlying structure of truth 
values, but we show that there are fuzzy automata such that there is a simulation or bisimulation of a given 
type between them, and there is not any crisp simulation/bisimulation of this type (cf. Example [Oj. 

The structure of the paper is as follows. In Section|2]we give definitions of basic notions and notation con- 
cerning fuzzy sets and relations, and in Section |3| we define basic notions and notation concerning fuzzy 
automata. In Section |4| we recall definitions of two types of simulations and four types of bisimulations 
introduced in lUsIl , and for any of them we define the corresponding function on the lattice of fuzzy relations. 
Section|5]contains our main results on the computing of the greatest simulations and bisimulations between 
fuzzy automata. Examples presented in Section|6|demonstrate the application of our algorithms and clarify 
relationships between different types of simulations and bisimulations. 
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It is worth noting that algorithms for deciding whether there are simulations and bisunulations between 
nondeterministic automata, and the computing of the greatest simulations and bisimulations, if they exist, 
were given in [14, 38]. Various algorithms for the computing of the greatest bisimulation equivalences on 
labelled transition systems can be found in ||2^|23,l36, 45, 55, 58]. 



2. Preliminaries 

The terminology and basic notions in this section are according to ll6l,[7l [T9ll20ll37ll . 

We will use complete residuated lattices as the structures of membership (truth) values. Residua ted 
lattices are a very general algebraic structure and generalize many algebras with very important applications 
(see for example |3. Bill' 111])- A residuated lattice is an algebra X = {L, A, V, ®, — >, 0, 1) such that 

(LI) (L, A, V, 0, 1) is a lattice with the least element and the greatest element 1, 
(L2) (L, ®, 1) is a commutative monoid with the unit 1, 

(L3) <8) and form an adjoint pair, i.e., they satisfy the adjunction property: for all x,y,z e L, 

x®y^z^x^y^z. (1) 

If, in addition, (L, A, V, 0, 1) is a complete lattice, then X is called a complete residuated lattice. Emphasizing 
their monoidal structure, in some sources residuated lattices are called integral, commutative, residuated 
/'-monoids [29]. 

The operations ® (called multiplication) and — > (called residuum) are intended for modeling the con- 
junction and implication of the corresponding logical calculus, and supremum (V) and infimum (/\) are 
intended for modeling of the existential and general quantifier, respectively. An operation <-> defined by 

X ^ y ^ {x ^ y) A{y ^ x), (2) 

called biresiduum (or biimplication), is used for modeling the equivalence of truth values. 

The most studied and applied structures of truth values, defined on the real unit interval [0, 1] with 
X A y = min(x, y) and x V y = max(x, y), are the Lukasiewicz structure (where x (8) y = max(x + y - 1, 0), 
X y = min(l - x + y, 1)), the Goguen (product) structure (x(g)y = x-y, x^y = lifx<y, and = y/x 
otherwise), and the Godel structure {x®y = min(x, y), x — > y = 1 if x < y, and = y otherwise). More generally, 
an algebra ([0, 1], A, V, ®, 0, 1) is a complete residuated lattice if and only if (g) is a left-continuous t-norm 
and the residuum is defined by x ^ y = Vl" ^ [0, 1] | m ® x < y\ (cf. |0])- Another important set of truth 
values is the set {flo, fli, . . . , fl„}, = aq < ■ ■ • < = 1, with a^^ai = flmax(*:+/-H,o) and a^ ^ a\ = flmin(n-S:+;,H)- 
A special case of the latter algebras is the two-element Boolean algebra of classical logic with the support 
(0,1). The only adjoint pair on the two-element Boolean algebra consists of the classical conjunction and 
implication operations. This structure of truth values we call the Boolean structure. A residuated lattice £. 
satisfying x^y-xAyis called a Heyting algebra, whereas a Heyting algebra satisfying the prelinearity 
axiom (x ^ y) V (y ^ x) = 1 is called a Godel algebra. If any finitelly generated subalgebra of a residuated 
lattice X is finite, then X is called locally finite. For example, every Godel algebra, and hence, the Godel 
structure, is locally finite, whereas the product structure is not locally finite. 

If X is a complete residuated lattice, then for all x, y, z e L and any {y,),e7 c L the following holds: 

X < y implies x ® z < y ® z, (3) 
X < y « X ^ y = 1, (4) 

X® Yy, = Y(x®y,), (5) 

!£j /el 

x®/\y, </\(x®y,), (6) 

For other properties of complete residuated lattices we refer to 
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In the further text £, will be a complete residuated lattice. A fuzzy subset of a set A over £, or simply a 
fuzzy subset of A, is any function from A into L. Ordinary crisp subsets of A are considered as fuzzy subsets 
of A taking membership values in the set {0, 1) c L. Let / and g be two fuzzy subsets of A. The equality of 
/ and g is defined as the usual equality of fxmctions, i.e., f = g ii and only if f{x) = g{x), for every x e A. 
The inclusion f ^ gis also defined pointwise: f ^ gif and only if f{x) ^ g{x), for every x g A. Endowed with 
this partial order the set ;F(A) of all fuzzy subsets of A forms a complete residuated lattice, in which the 
meet (intersection) Aiei/i ^r^d the join (union) V/ei.^' of an arbitrary family {fjiei of fuzzy subsets of A are 
functions from A into L defined by 



and the product / (gi ^ is a fuzzy subset defined by / (gi g{x) = f{x) (8 g{x), for every x e A. 

Let A and B be non-empty sets. A fuzzy relation between sets A and B (or & fuzzy relation from A to B) is any 
function from A X B into L, that is to say, any fuzzy subset of A X B, and the equality, inclusion (ordering), 
joins and meets of fuzzy relations are defined as for fuzzy sets. In particular, a fuzzy relation on a set A is any 
function from AxA into L, i.e., any fuzzy subset of Ax A. To highlight the difference between fuzzy relations 
between two sets and those on a set, fuzzy relations between two sets will be usually denoted by small 
Greek letters, and fuzzy relations on a set by capital Latin letters. The set of all fuzzy relations from A to B 
will be denoted by "RiA, B), and the set of all fuzzy relations on a set A will be denoted by '/?(A). The converse 
(in some sources called inverse or transpose) of a fuzzy relation (p e ??(A, B) is a fuzzy relation (p~^ e 'R{B, A) 
defined by (p~^{b,a) = (p{a, b), for all a e A and b e B. A crisp relation is a fuzzy relation which takes values 
only in the set {0, 1|, and if ^ is a crisp relation of A to B, then expressions "(p{a, b) = 1" and "(a, b) e cp" will 
have the same meaning. 

For non-empty sets A, B and C, and fuzzy relations 93 e ??(A, B) and ip e 'R{B, C), their composition cp oip 
is an fuzzy relation from ??(A, C) defined by 

{(p o xp){a, c) = Y (pia, b) ® \p{b, c), (7) 

for all a 6 A and c e C. If cp and \p are crisp relations, then ^ o is an ordinary composition of relations, i.e., 

(poTl> = {{a, c) e A x C I (3fc 6 B) (fl, fc) e 9? & {b, c) e i/^} , 

and if (p and i/' are functions, then ip o is an ordinary composition of functions, i.e., {(p o \p){a) = ip{(p{a)), for 
every a e A. Next, if / e !F(A), <p e '/?(A, B) and ^ e ^(B), the compositions f °(p and <p o ^ are fuzzy subsets 
of B and A, respectively, which are defined by 

(/ ° (pm = V /(«) ® <p{a, b), {cp o ^)(fl) = Y (p(fl, [;) ® g{b), (8) 

«€A ij€B 

for every a e A and b e B. 

In particular, for f,g& !^(A) we write 

fog = \J f{a)®g{a). (9) 

The value fog can be interpreted as the "degree of overlapping" of / and g. In particular, if / and g are 
crisp sets and ip is a crisp relation, then 

f o (p = [b e B \ {3a e f) {a,b) e (p] , cp o g = [a & A \ {3b & g) {a,b) e cp] . 

Let A, B, C and D be non-empty sets. Then for any (p\ e ??(A, B), e ??(B, C) and (ps e ??(C, D) we have 



{(pi o 932) o 933 = o (9J2 o 933), 
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(10) 



and for (po e ??(A, B), (pi, (p2 e 'R(B, C) and (pj e ??(C, D) we have that 



(j9i ^ (j92 implies (p^ ^ 'P2 > (Po ° (pi ^ (po ° 'Pi/ and (p\ o cp^ ^ o (pj. 



(11) 



Further, for any (p e 'R(A, B), i/^ e 'R(B, C), / e F(A), ^ e T(B) and e ;r(C) we can easily verify that 



We note that if A, B and C are finite sets of cardinality |A| - k, |B| = m and |C| = n, then ip e "^(A, B) and 
1/^ e "^(B, C) can be treated as kxm and m x n fuzzy matrices over X, and cp o 1^ is the matrix product. Analo- 
gously, for / e IF (A) and g e F(B) we can treat f o cp as the product of a 1 x A: matrix / and akxm matrix cp 
(vector-matrix product), (p o g as the product of an A: X m matrix (p and an m X 1 matrix g' , the transpose of g 
(matrix-vector product), and f o g as the scalar product of vectors / and g. 



3. Fuzzy automata 

In this paper we study fuzzy automata with membership values in complete residuated lattices. 

In the rest of the paper, if not noted otherwise, let X be a complete residuated lattice and let X be an (finite) 
alphabet. A fuzzy automaton over X. and X , or simply a fuzzy automaton, is a quadruple J?l = (A, 6^, a^, t^), 
where A is a non-empty set, called the set of states, 6^ :AxXxA—>L is a fuzzy subset of A x X x A, called 
the fuzzy transition function, and : A ^ L and t"^ : A ^ L are the fuzzy subsets of A, called the fuzzy set of 
initial states and the fuzzy set of terminal states, respectively. We can interpret 6^ {a, x, b) as the degree to which 
an input letter x e X causes a transition from a state a e A into a state b e A, whereas we can interpret o^{a) 
and T^{a) as the degrees to which a is respectively an input state and a terminal state. 

Many results presented in this paper can also be proved for fuzzy automata with infinite sets of states. 
However, here we deal with effective algorithms for computing the greatest simulations and bisimulations, 
and for that reason we will assume that the sets of states of the considered fuzzy automata are always/zVizfe. 

Let X* denote the free monoid over the alphabet X, and let e e X* be the empty word. The fimction 6"^ 
can be extended up to a function 6^ : A x X* x A ^ L as follows: It a,b e A, then 



{f o (p) o ip - f o [(p o ip), {f o (p) o g - f o {(p o g), {(p o ip) o h = (p o (ip o h) 



(12) 



and consequently, the parentheses in |[T2t can be omitted, as well as the parentheses in JTOll . 
Finally, for all (p, (pi e '/?(A, B) (z e I) and ip, ipi e 'R(B, C) (z e I) we have that 

{(p o \p)~^ = \p~^ o (p~^. 



(13) 
(14) 




(15) 




(16) 



and if a,b e A, u e X* and x eX, then 




(17) 



By Q and Theorem 3.1 ^ (see also Isl-lsl), we have that 



6^ {a, uv, b) = \J 6^ {a, u, c) (g) 6f (c, v, b). 
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(18) 



for all a,b e A and ii, v e X*, i.e., if if = xj ■ • ■ Xn, for xi,...,x„ eX, then 

6^{a,iv,b)= \J 6^(fl,xi,ci)®6^(ci,X2,C2)® •••®6'^(c„_i,x„,J;). (19) 

(ci,...,c„_i)€A"-i 

Intuitively, the product 6'^{a,xi,ci) (gi d^{ci,X2,C2) (8i • • • ® 6^{c„-i,x„,b) represents the degree to which the 
input word w causes a transition from a state a into a state through the sequence of intermediate states 
Ci, . . . , c„_i e A, and w, b) represents the supremum of degrees of all possible transitions from a into 
b caused by iv. Also, we can visualize a fuzzy finite automaton J?l representing it as a labelled directed 
graph whose nodes are states of and an edge from a node a into a node b is labelled by pairs of the form 
x/5^{a, X, b), for any x e X. 

It is worth noting that fuzzy automata are a natural generalization of nondeterministic and deterministic 
automata. Namely, if is a crisp subset of A X X X A, that is, 6^ : A X X X A {0, 1 }, and and t"^ are crisp 
subsets of A, then J[ is an ordinary nondeterministic automaton. In other words, nondeterministic automata 
are fuzzy automata over the Boolean structure. They will also be called Boolean automata. If is a function 
of A X X into A, is a one-element crisp subset of A, that is, - {floL for some aq ^ ^/ arid is a fuzzy 
subset of A, then J?I is called a deterministic fuzzy automaton, and it is denoted by ^ = (A, 6'*,flo/T'*)- hi 
iislll^ the name crisp-deterministic was used. For more information on deterministic fuzzy automata we 



refer to |a,|30H31|,|Ml3S,^]- Evidently, if b^ is a crisp subset of A x X x A, or a function of A x X into A, then 
6^ is also a crisp subset of A x X* x A, or a function of A x X* into A, respectively. A deterministic fuzzy 
automaton JTI = (A, 6'*, aq/ t^)/ where t"^ is a crisp subset of A, is an ordinary deterministic automaton. 
If for any m e X* we define a fuzzy relation 6^ on A by 

6^(fl,b) = 6^'(fl,u,&), (20) 
for all a,b e A, called the fuzzy transition relation determined by u, then dTSl l can be written as 

bt = b^o5t (21) 
for all u,v e X*. 

The reverse fuzzy automaton of a fuzzy automaton J?l = (A, 5"^, a^, t^) is defined as the fuzzy automaton 
^ = (A, 6'^, cj'^, T^) whose fuzzy transition function and fuzzy sets of initial and terminal states are defined 
by B^{ai,x,a2) = b^{a2,x,ai) for all fli,fl2 £ A and x & X, a'^ = and = a^. In other words, = (6^^)"^, 
for each x e X. 

A fuzzy language in X* over X, or briefly a fuzzy language, is any fuzzy subset of X*, i.e., any function 
from X* into L. A fuzzy language recognized by a fuzzy automaton J{ = {A,b^ ,a^,T^), denoted as L{J{), is a 
fuzzy language in !F(X*) defined by 

L{^){u) = Y a^(fl) ® bt{a, u, b) ® T^{b), (22) 

a,b€A 

or equivalently, 

L(:?l)(e) = a^oT^ 

L(^)(m) = o o o . . . o o T^, 

for any u = X1X2 • • . x„ e X^, where xj, X2/ ■ • ■ / £ X. In other words, the equality l|22|l means that the mem- 
bership degree of the word u to the fuzzy language L{J{) is equal to the degree to which J{ recognizes or 
accepts the word u. Using notation from (|8}, and the second equality in ((T2I 1, we can state (|22] l as 

L{Ji){u) = o 6^ o T^. (24) 

Fuzzy automata J{ and S are called language-equivalent, or sometimes just equivalent, if L(J?l) = i^(S). 

For more information on fuzzy automata with membership values in complete residuated lattices we 
refer to the recent papers lll6, 17, 3lll33,l33,l52l-l53 . l63.l66l - [70^ and for information on fuzzy automata over 



related types of lattices we refer to The most complete overview of the results concerning the 

"classical" fuzzy automata taking membership values in the Godel or the Goguen structure can be found 
in the book JH. 
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4. Simulations and bisimulations between fuzzy automata 

First we recall definitions of two types of simulations and four types of bisimulations introduced in f\Bi\ . 
Let Jl = (A,6^,0^,T^) and S = {B,6'^,a'^,T'^) be fuzzy automata, and let cp e 'R{A,B) be a non-empty 
fuzzy relation. We call <p a forward simulation if it satisfies 

< o (p-\ (fs-1 
(p~^ 06^ ^ 6^ o (p~^, for every xeX, {fs-1 
(p-^ o < (/s-3 

and a backward simulation if 



A 



6^ o cp ^ (p o 5^, for every x eX, (bs-2 
ocp^a^. (bs-3 

Furthermore, we call cp a forward bisimulation if both <p and ^"^ are forward simulations, i.e., if cp satisfies 

a^<a^o (p-\ < a"^ o cp, {fb-1 

(p~^ o 6^ < 6^ o (p-^, (pobl^b^ o (p, for every xeX, (fb-1 

and a backward bisimulation, if both ^ and <p~^ are backward simulations, i.e., if <p satisfies 

T^^ < ip o T^, < ^"^ o T^, {bb-1 

o(p ^(podl, bl o < (p-1 o 6^, for every xeX, {bb-1 
a^o(p^ a^, o cp-i < o^. {bb-3 

Also, if cp is a forward simulation and (p~^ is a backward simulation, i.e., if cp satisfies 

^""^ 06^ = ° 'P~^ for every x e X, if bb-1 

then cp is called a forward-backward bisimulation, and if (p is a backward simulation and (p~^ is a forward 
simulation, i.e., if 

cr^ < o ^, t"^ < ^ o T^, (^'/^'-i: 

b^ o (p = (p o b^, for every x e X, (bfb-l 

a^ocp^a^ (poT^^T^. {bfb-3 

then (p is called a backward-forward bisimulation. 

For the sake of simplicity, we will call cp just a simulation if cp is either a forward or a backward simulation, 
and just a bisimulation if cp is any of the four types of bisimulations defined above. Moreover, forward 
and backward bisimulations will be called homotypic, whereas backward-forward and forward-backward 
bisimulations will be called heterotypic. 

The meaning of forward and backward simulations can be best explained in the case when ^ and S are 
nondeterministic (Boolean) automata. For this purpose we will use the diagram shown in Figure 1 . Let (p be a 
forward simulation between ^ and S and let ao,ai, . . ., a„ be an arbitrary successful run of the automaton ^ 
on a word u = X\X2 • • • x„ {x\, X2,...,x„ e X), i.e., a sequence of states of M such that aq ^ {^k, fl/t+i) e ^x^+i' 
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Figiire 1: Forward and backward simulation 



forO ^ < n-l,andfl„ e . According to ( fs-l\ , there is an initial state bp e such that (flo/^o) £ (p. Suppose 



that for some A;, < A: < n - 1, we have built a sequence of states bo, bi, . . .,bjf such that b,) e 6f and 
(fl„b,) e (p, for each;, 1 < f ^ fc. Then {bk,ak+i) e ip"^ o 6^,^^^ , and by l /s-2"| we obtain that {bk,ak+i) e 6^^,^^ °(p~^, 



so there exists fc^+i £ B such that (fc/^, b/c+j) e Sj:^^^^ and e cp. Therefore, we have successively built 

a sequence bQ,bi, ...,b„ of states of B such that bo e "J^/ {bk,bk+i) e S^^^ for every fc, < fc < n - 1, and 
(uk, bk) e <p, for each A:, < fc < n. Moreover, by ( /s-3| we obtain that b„ e t^. Thus, the sequence bo, bi,...,b„ 



is a successful run of the automaton S on the word ii which simulates the original run ao,ai, ■ . . ,a„ oi J{ 
on u. In contrast to forward simulations, where we build the sequence bQ, b\, . . ., b„ moving forward, starting 
with bp and ending with bn, in the case of backward simulations we build this sequence moving backward, 
starting with b,, and ending with I^q. In a similar way we can understand forward and backward simulations 
between arbitrary fuzzy automata, taking into account degrees of possibility of transitions and degrees of 
relationship. 

In numerous papers dealing with simulations and bisimulations mostly forward simulations and for- 
ward bisimulations have been studied. They have been usually called just simulations and bisimulations, 
or strong simulations and strong bisimulations (cf. ll42l,l43l.l56l] V and the greatest bisimulation equivalence has 
been usually called a bisimilarity . Distinction between forward and backward simulations, and forward and 
backward bisimulations, has been made, for instance, in ESEiliO] (for various kinds of automata), but less 



or more these concepts differ from the concepts having the same name which are considered here. More sim- 
ilar to our concepts of forward and backward simulations and bisimulations are those studied in [9], and in 
iIztIEsIi (for tree automata). Moreover, backward-forward bisimulations have been discussed in the context 
of weighted automata in |lfl|i|l3,El|2i]. 

Given fuzzy automata J?I = (A, b^, a^, t^) and B = (B, 5^, , t^), a fuzzy relation (p e 'K(A, B) is a back- 
ward simulation between fuzzy automata J\ and S if and only if it is a forward simulation between the 
reverse fuzzy automata J{ and S. Therefore, forward and backward simulations, forward and backward 
bisimulations, and backward-forward and forward-backward bisimulations, are mutually dual concepts, 
what means that for any statement on some of these concepts which is universally valid (valid for all fuzzy 
automata) there is the corresponding universally valid statement on its dual concept. 

Now we introduce several other notions and notation that will be used in future work. 

For non-empty sets A and B and fuzzy subsets f] e T^iA) and E, e 7^(B), fuzzy relations j] ^ e 'R{A, B) 
and f] <— i£ e '^(A, B) are defined as follows 

{^^^){a,b) = {iia) ^ m), (25) 
(r] ^ .£)(«,&) = (.£(&) ^ ,](«)), (26) 

for arbitrary a e A and b e B. Let us note that ?] <— ^ = (<^ — > 7])~^. 
We have the following. 
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Lemma 4.1. Let A and B be non-empty sets and let rj e IF (A) and £, e !F(B). 

(a) The set of all solutions to the inequality rjo x ^ where x is an unknown fuzzy relation between A and B, is 
the principal ideal of'R{A, B) generated by the fuzzy relation rj 

(b) The set of all solutions to the inequality x° ^ ?]/ where x is an unknown fuzzy relation between A and B, is 
the principal ideal of'R{A, B) generated by the fuzzy relation ri *— 

Proof. These are the well-known results by E. Sanchez (cf . llsil - 61 1)- □ 

It is worth noting that (?] £,) A (;] <— <5) = ?] <5, where rj <-> (J is a fuzzy relation between A and B 
defined by 

{!] ^ Q{a,b) ^ {i]{a) ^ mi (27) 

for arbitrary a e A and b e B. 

Next, let A and B be non-empty sets and let a e !/?(A), jS e 'R{B) and (p e B). The right residual of <p by 
a is a fuzzy relation (p/a e 'R(A, B) defined by 

{(p/a){a, b)= /\i a{a', a) ^ <p(fl', b) ), (28) 

for all fl e A and b eB, and the left residual of ^ by /3 is a fuzzy relation ^\jS e ^{A, B) defined by 

{(p\^){a, b)= /\( (iib, b') ^ <p(a, V) ), (29) 

for all fl e A and b eB.ln the case when A = B, these two concepts become the well-known concepts of right 
and left residuals of fuzzy relations on a set (cf. [33]). 
We also have the following. 

Lemma 4.2. Let A and B be non-empty sets and let a e 'R(A), jS e '??(B) and cp e KiA, B). 

(a) The set of all solutions to the inequality ao x ^(p, where x is an unknown fuzzy relation between A and B, is 
the principal ideal of'R{A, B) generated by the right residual (p/a ofofcp by a. 

(b) The set of all solutions to the inequality X ° § ^ 'P^ where x is an unknown fuzzy relation between A and B, is 
the principal ideal of'R{A, B) generated by the left residual (p\fi ofofcp by f}. 



Proof. These are also results by E. Sanchez (cf. Il59l - l61ll ). □ 

As we said in the introduction, the problem of deciding whether there is a simulation or bisimulation of 
a given type between fuzzy automata, and the problem of computing the greatest simulation/bisimulation 
of this type, will be reduced to the problem of computing the greatest post-fixed point, contained in a given 
fuzzy relation, of an appropriate isotone function on the lattice of fuzzy relations. For this purpose we 
define the following fuzzy relations and functions on the lattice of fuzzy relations. 

Let ^ = (A, 5^, o^, T^) and B = (B, 6^, a^, t^) be fuzzy automata. We define fuzzy relations f" e !/?(A, B), 
for w e j/s, bs, fb, bb, fbb, bfb], in the following way: 

i/,/" = ^ T^, (30) 

= ff^ ^ (jfi, (31) 

^fl = (T^ ^ T^) A (T^ ^ T^) = « T^ (32) 
^bb ^ (^A ^ ^B^ ^ (^A ^ ^B) ^^A^ ^B^ (33) 

^fbb = (^A ^ ^B) ^ (^A ^ ^B)^ (34) 

^bfb ^ (^A ^ ^B) ^ (^A ^ ^By (35) 



Moreover, we define functions cf)''" : '??(A, B) "RiA, B), for w e {fs, bs,fb, bb,fbb, bfb], as follows: 

cpf%a) = /\moa-')\5i]-\ (36) 

.vex 

cp''%a) = /\{ao6^,)/dt (37) 

vex 

(p^%) = /\[{l>x ° «"')\6v A [(6^ o a)\6^] = (pf^a) A [cpf%a-')]-\ (38) 

.vex 

cp'Hcc) = /\[{a o 6^)/6^] A [(«-! o 6^)/6f ]-i = cp^'^ia) A [c/)^X«"')]"\ (39) 

.vex 

(pf'\a) = /\[(6f o a-i)\6^]-i A [(a'^ o S^)/^^]-! = (|,/-^(a) A [(/)''^(a-i)]-i, (40) 

vex 

<p'^\a) = /\[{ao 6f )/6^] a [(5^ o a)\6f] = c/)^X«) A [cpHa-')]-\ (41) 

vex 

for any a e '??(A,B). Notice that in the expression "(p^''{a~^)" {lu e {fs, bs]) we denote by (p^ a function from 
'R(B,A) into itself. 

The next theorem provides equivalent forms of the second and third conditions in the definitions of simu- 
lations and bisimulations. 

Theorem 4.3. Let J{ = {A, 6^, a^, t^) and S = (B, 6^, a^, t^) be fuzzy automata and w e {fs, bs, fb, bb, fbb, bfb). 
A fuzzy relation (p e 'R{A,B) satisfies conditions {w-2) and {w-3) if and only if it satisfies 

cp < f"{(p), (p < f. (42) 

Proof. We will prove only the case iv - fs. The assertion concerning the case zv = bs follows by the duality, 
and according to equations ll32t - ll35] l and ll38l l- (|4T] |, all other assertions can be obtained by the first two. 

Consider an arbitrary (p e ^^{A, B). According to Lemma [4.1f b), (p satisfies condition (/s-3) if and only if 
(p~^ < <— T"^ = (t"^ T^)"^, which is equivalent to cp < t"^ ^ = \pf^. Therefore, cp satisfies (/s-3) if and 
only if (p < \pf^. 

On the other hand, (p satisfies (/s-2) if and only if 

(p-Hb,a)<^6^(a,a') < {d'^, o (p-^)(b,a'), 
for all a,a' e A, b £ B and x e X. According to the adjunction property, this is equivalent to 
cp-Hb,a) ^ l\[d^{a,a') ^ {6^ o <p-'ib,a'))] = {{5^ o cp-')\5^){b,a) 

a'€A 

for all fl e A, e B and x e X, which is further equivalent to 
cp{a,b) < /\[{6^ o cp-')\d^]-Ha,b) = {cpH<p))(a,b) 

vex 

for all fl e A and b e B. Therefore, cp satisfies (/s-2) if and only if ip < (pf^{(p). 

Now we conclude that a fuzzy relation cp e 'R{A, B) satisfies (/s-2) and (/s-3) if and only if it satisfies 
(for w - fs), which was to be proved. □ 

5. Computation of the greatest simulations and bisimulations 

We will provide a method for computing the greatest simulations and bisimulations between fuzzy auto- 
mata adapting the method developed in [33] for computing the greatest post-fixed points of an isotone 
function on the lattice of fuzzy relations on a single set. 
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Let A and B be non-empty sets and let : ??(A, B) 'R(A, B) be an isotone function, i.e., let a < jS implies 
<\){ci) < for all a, jS e '??(A, B). A fuzzy relation a e ^^(A, B) is called a post-fixed point of c/) if a < (^(a). The 
well-known Knaster-Tarski fixed point theorem (stated and proved in a more general context, for complete 
lattices) asserts that the set of all post-fixed points of (p form a complete lattice (cf. |57]). Moreover, for any 
fuzzy relation ijj e 'R(A, B) we have that the set of all post-fixed points of (p contained in ijj is also a complete 
lattice. According to Theorem l4.31 our main task is to find an effective procedure for computing the greatest 
post-fixed point of the function cp'" contained in the fuzzy relation jp^, for each iv e {/s, bs,fb, bb,fbf, bfb]. 

It should be noted that the set of all post-fixed points of an isotone function on a complete lattice is always 
non-empty, because it contains the least element of this complete lattice. However, this set may consist only 
of that single element. In our case, when we are working with a lattice of fuzzy relations, the empty relation 
may be the only post-fixed point, whereas we have defined simulations and bisimulations to be non-empty 
fuzzy relations. This requirement is necessary because the empty relation can not satisfy the condition (iv-l), 
unless the fuzzy set of initial states or the fuzzy set of terminal states is also empty. Therefore, our task is 
actually to find an effective procedure for deciding whether there is a non-empty post-fixed point of (p^' 
contained in ip'", and if it exists, then find the greatest one. 

Let (p : '??(A, B) —> 'R(A, B) be an isotone function and ip e ??(A, B). We define a sequence {(pk]k€N of fuzzy 
relations from 'R(A, B) by 

cpi = ip, cpk+i = (pk ^ (p{(Pk), for each A: e N. (43) 

The sequence {<pk]k€M is obviously descending. If we denote by <p the greatest post-fixed point of (p contained 
in \p, we can easily verify that 



Now two very important questions arise. First, under what conditions the equality holds in (|44t ? Even 
more important question is: under what conditions the sequence {cptliceN is finite? If this sequence is finite, 
then it is not hard to show that there exists A: e N such that (pk = (p,,,, for every m ^ k, i.e., there exists e N 
such that the sequence stabilizes on (p^. We can recognize that the sequence has stabilized when we find the 
smallest fc e N such that qj^ - fk+i- Iri this case <p = (pk, and we have an algorithm which computes ^ in a 
finite number of steps. 

Some conditions under which equality holds in (|44] l or the sequence is finite were found in [33], for 
fuzzy relations on a single set. It is easy to show that the same results are also valid for fuzzy relations 
between two sets. In the sequel we present these results. 

First we note that a sequence {(pk]k€N of fuzzy relations from 'R(A, B) is finite if and only it it is image- 
finite, which means that the set U/ceN I^l?**:) is finite. Next, the function <p : 'R(A,B) 'R(A,B) is called 
image-localized if there exists a finite X c L such that for each fuzzy relation (p e '/?(A, B) we have 



where (K U Im(^)) denotes the subalgebra of £. generated by the set K U Im(^). Such K will be called a 
localization set of the function (p. 

Theorem analogous to the following theorem was proved in \ for fuzzy relations on a single set, but its 
proof can be easily transformed into the proof of the corresponding theorem concerning fuzzy relations 
between two sets. 

Theorem 5.1. Let the function (p be image-localized, let K be its localization set, let Tp e ??(A, B), and let {(pk]k€f^ be 
a sequence of fuzzy relations in 'R{A,B) defined by i43i . Then 




(44) 



lm{(p{q))) c(KU Im((p)), 



(45) 




(46) 



If, moreover, (K U lm{\p)) is a finite subalgebra of then the sequence {(pk]k€n is finite. 
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Going back now to the functions (/)'"', for zv e {fs, bs, fb, bb, fbb, bfb], we prove the following. 

Theorem 5.2. Let J{ = {A, 6^, a^, t^) and S = (B, 6^, a^, t^) be arbitrary automata. 
For any w e [fs, bs,fb, bb,fbb, bfb] the function (p^" is isotone and image-localized. 

Proof. We will prove only the case lu - fs. The assertion concerning the case w = bs follows by the duality, 
and according to equations ll32t-ll35]l and ll38ll- (|4T] |, all other assertions can be derived from the first two. 

Let <pi,(p2 G fi(-A,B) be fuzzy relations such that (pi < q)2, and consider the following systems of fuzzy 
relation inequalities: 

o 5^ < 6^ o (p-\ X e X; (47) 
X"^ o 6^ < 6f o (p-\ X e X. (48) 

where x ^ fii^, B) is an unknown fuzzy relation. Using Lemma [4.2l (b) and the definition of an inverse rela- 
tion, it can be easily shown that the set of all solutions to system ll47t (resp. (|48] |) form a principal ideal of 
'R(A,B) generated by (|)-^'''(^i) (resp. c/)^'X(p2))- Since for each x e Xwehavethat6^o^~i < 6^o(p~i, we conclude 
that every solution to (|47ll is a solution to (|48] l. Consequently, (pf^{(pi) is a solution to (|48l l, so (pf^cpi) < (pf'^{(p2). 
Therefore, we have proved that (p^'^ is an isotone function. 

Next, let K = Uxexi^^i^x) Iin(^x)) arid let (p e 'R{A, B) be an arbitrary fuzzy relation. It is evident that 
lin((pf^{(p)) c (X U lm(<p)), and since the input alphabet X is finite, then K is also finite. This confirms that 
the function (pf^ is image-localized. □ 

Now we are ready to prove the main result of this paper, which provides algorithms for deciding whether 
there is a simulation or bisimulation of a given type between fuzzy automata, and for computing the greatest 
simulations and bisimulations, when they exist. 

Theorem 5.3. Let J{ = (A,6^,a^,T^) and S = (B,6S,a^,T^) be fuzzy automata, let w e {fs,bs, fb,bb, fbb,bfb], 
and let a sequence {(pk]k&i of fuzzy relations from ^{A, B) be defined by 

(pi = ip'", (pk+i = (pk^ (p^''i<Pk), for each ke'N. (49) 

If{lm{}p^") U U.v6x( ) U ^^i^x) )) a finite subalgebra of then the following is true: 

(a) the sequence [(pk\k€K is finite and descending, and there is the least natural number k such that (p^ = (pk+ii 

(b) (pk is the greatest fuzzy relation in ??(A, B) which satisfies {w-2) and {w-3); 

(c) if(pk satisfies {w-\), then it is the greatest fuzzy relation in 'R{A,B) which satisfies {w-\), {w-2) and {w-3); 

(d) if(pk does not satisfy {w-1), then there is no any fuzzy relation in 'R{A, B) which satisfies {w-1), {w-2) and {w-3). 

Proof. We will prove only the case w = fs. All other cases can be proved in a similar manner. 
So, let {lm{^>f') U Ux€x( Im(6^) U Im(6^) )) be a finite subalgebra of X. 

(a) According to Theorems l5.2l and l5.1[ the sequence {(pk]k(iN is finite and descending, so there are fc, m e IN 
such that (pk = (pk+iru whence cpk+i ^ (pk ~ (Pk+m ^ (Pk+i- Thus, there is A: e N such that (pk - (pk+i, and conse- 
quently, there is the least natural number having this property. 

(b) By (Pit = <Pk+i = (Pk^ 'p^'^i'Pk) we obtain that (pk < (p^^{(pk), and also, <pk ^ <pi = ^f'^. Therefore, by Theo- 
rem l4.3l it follows that (pk satisfies (/s-2) and (/s-3). 

Let (p e ^{A, B) be an arbitrary fuzzy relation which satisfies (/s-2) and (/s-3). As we have already noted, 
(p satisfies (/s-3) if and only if < \pf^ = (pi. Next, suppose that (p < (p„, for some « e N. Then for every 
X e X we have that (p~^ ° ^ ^.t ° ^ ° 'Pn^' ^rid according to Lemma l42l (b), (p~^ < (6f o (p~i) \ 5^, 
i.e., (p ^ [{5x o (p~i) \ 6^]"^ = (pf^{(p„). Therefore, <p ^<p„ A (pf%(p„) = (p„+i. Now, by induction we obtain that 
(p < (p„, for every n e N, and hence, (p < qjk- This means that (pk is the greatest fuzzy relation in 'R{A, B) 
satisfying (/s-2) and (/s-3). 

(c) This follows immediately from (b). 
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(d) Suppose that (pj^ does not satisfy (fs-l). Let cp e '/?(A, B) be an arbitrary fuzzy relation which satisfies 
(/s-1), (/s-2) and (/s-3). According to (b) of this theorem, cp < (pi^, so we have that < a^o<p~'^ < °'Pk^ ■ ^^t, 
this contradicts our starting assumption that (p^ does not satisfy (/s-1). Hence, we conclude that there is no 
any fuzzy relation in 'R(A, B) which satisfies (/s-1), (/s-2) and (/s-3). □ 

Next, we will consider the case when X - (L, A, V, ®, 0, 1) is a complete residuated lattice satisfying 
the following conditions: 

xv{/\y,) = /\{xvy,l (50) 

i€l iel 

x®(/\y,) = /\(x®y/), (51) 



/€/ iel 



for all X e L and {i/,j,ez £ L. Let us note that if X = ([0, 1], A, V, ®, 0, 1), where [0, 1] is the real unit interval 
and ® is a left-continuous t-norm on [0, 1], then (ISOl l follows immediately by linearity of X, and X satisfies 
JSTIl if and only if ® is a continuous t-norm, i.e., if and only if X is a BL-algebra (cf . @0])- Therefore, conditions 
iH and ^ hold for every BL-algebra on the real unit interval. In particular, the Lukasiewicz, Goguen 
(product) and Godel structures fulfill l|50|l and (Sit . 
Under these conditions we have the following. 

Theorem 5.4. Let J{ = (A, 6^, ct^, t^) and S = (B, 6^, , t^) be fuzzy automata, let w e {/s, bs, fb, bb, fbb, bfb], 
let {(pk]k€K l^s ^ sequence of fuzzy relations from 'R{A, B) defined by f49[ , and let 

(p = (52) 

If £. is a complete residuated lattice satisfying (|50t and | |5I] |, then the following is true: 

(a) (p is the greatest fuzzy relation in 'R{A,B) which satisfies {w-2) and (if-3); 

(b) if(p satisfies {w-1), then it is the greatest fuzzy relation in '??(A, B) which satisfies {w-1), {w-2) and {w-3); 

(c) if(p does not satisfy {w-1), then there is no any fuzzy relation in ^{A, B) which satisfies {w-1), {w-2) and {w-3). 

Proof. We will prove only the case w = fs. All other cases can be proved similarly. 

So, let X be a complete residuated lattice satisfying llSOt and dSTT l. First, notice that it has been shown in 
iIitIi that if (l50l holds, then for all descending sequences [xk]k€Hf Ij/fclteN £ i we have 

/\{xkyyk) = (/\xk)y(/\yk)- (53) 

k€M k<m k€M 

(a) For arbitrary x e X, a e A and b eB we have that 



\{6lo(pl^)Yb,a)= \{dlo(pl^){b,a) = \\AJ 

teN ' teN kaM^b'aB ' 

= \j{/\dl{b,b')®cp-\b',a)\ (by®) 

= \j[bl{h,V)®il\cpl\b',a)^ (bydSB) 
bl{b, b') ® (p-\b', fl)) = {bl o (p-^){b, a), 



13 



which means that 



i-eN 

for every x e X. The use of condition (|53) is justified by the facts that B is finite, and that [<p'^^{b' ,a)]k€M is a 
descending sequence, so {6f (b, h') ® (p'^^{b' , a))te]N is also a descending sequence. 
Now, for all X e X and /c e N we have that 

which is equivalent to 

As the last inequality holds for every e N, we have that 

for every x e X. Therefore, cp satisfies (/s-2). Moreover, cp ^cpi = }pf^, so cp also satisfies (/s-3). 

Next, let a e '/?(A, B) be an arbitrary fuzzy relation satisfying (/s-2) and (/s-3). According to Theorem l4.31 
a < (p^'^{a) and a < i/;-'^^ = epi. By induction we can easily prove that a cpk for every k e N, and therefore, 
a K (p. This means that cp is the greatest fuzzy relation in 'R(A, B) which satisfies (/s-2) and (/s-3). 

The assertion (b) follows immediately from (a), whereas the assertion (c) can be proved in the same way 
as the assertion (d) of Theorem l5.3l □ 

In some situations we do not need simulations and bisimulations that are fuzzy relations, but those that 
are ordinary crisp relations. Moreover, in cases where our algorithms for computing the greatest simulations 
and bisimulations fail to terminate in a finite number of steps, we can search for the greatest crisp simulations 
and bisimulations. They can be understood as a kind of "approximations" of the greatest fuzzy simulations 
and bisimulations. Here we show that the above given algorithms for computing the greatest fuzzy simu- 
lations and bisimulations can be modified to compute the greatest crisp simulations and bisimulations. The 
new algorithms terminate in a finite number of steps, independently of the properties of the underlying 
structure of truth values. In the next section we will give an example which shows that the greatest crisp 
simulations and bisimulations can not be obtained simply by taking the crisp parts of the greatest fuzzy sim- 
ulations and bisimulations. In fact, our example shows that there may be a fuzzy simulation/bisimulation 
of a given type between two fuzzy automata, but there is not any crisp simulation/bisimulation of this type 
between them. 

Let A and B be non-empty finite sets, and let 'R.'^{A, B) denote the set of all crisp relations from 'R(A, B). It 
is not hard to verify that 'R'^{A, B) is a complete sublattice of ??(A, B), i.e., the meet and the join in '??(A, B) of an 
arbitrary family of crisp relations from '/?'^(A, B) are also crisp relations (in fact, they coincide with the ordi- 
nary intersection and union of crisp relations). Moreover, for each fuzzy relation <p e 'R{A, B) we have that 
cp'^ e 'R'^{A, B), where (p'^ denotes the crisp part of a fuzzy relation cp (in some sources called the kernel of cp), 
i.e., a function : A x B ^ {0, 1| defined by <p'^{a,b) = 1, if (p{a,h) = 1, and (p^{a,h) = 0, if (p(a,b) < 1, for 
arbitrary fl e A and b e B. Equivalently, (p'^ is considered as an ordinary crisp relation between A and B given 
by (p' = {(a, b)eAxB\ (p{a, h) = 1). 

For each fimction (p : 'R(A, B) ^ ??(A, B) we define a function cp'' : 'R%A, B) ^ 'R'iA, B) by 

(P'^lcp) = {(p{(p)Y, for any cp e 'R'iA, B). 

If (p is isotone, then it can be easily shown that (p'^ is also an isotone function. 
We have that the following is true. 

14 



Proposition 5.5. Let A and B be non-empty finite sets, let cp : ??(A, B) 'R{A,B) be an isotone function and let 
\p e ??(A, B) be a given fuzzy relation. A crisp relation q e 'R'^{A,B) is the greatest crisp solution in 'R{A,B) to the 
system 

X < (pixl X<4'' (54) 

if and only if it is the greatest solution in 'R'^{A, B) to the system 

^ < (P'{Q, ^ < 4>', (55) 

where x is an unknown fuzzy relation and E. is an unknown crisp relation. 
Furthermore, a sequence {p^ljteiN £ "J^iA, B) defined by 

Qi = ip", gk+1 = Qkfy (p'iek), for every fc e N, (56) 

is a finite descending sequence of crisp relations, and the least member of this sequence is the greatest solution to the 
system l|55| in 'R''{A, B). 

Proof. The proof of this proposition can be obtained simply by translating the proof of Theorem 5.8. fsS*] to 
the case of relations between the two sets. □ 

Taking (p to be any of the functions (p™ , for w e [fs, bs, fb, bb, fbb, bfb]. Prop osition l5 . 51 gives algorithms for 
deciding whether there is a crisp simulation or bisimulation of a given type between fuzzy automata, and for 
computing the greatest crisp simulations and bisimulations, when they exist. As we have seen in Proposi- 
tion l5.5l these algorithms always terminate in a finite number of steps, independently of the properties of the 
underlying structure of truth values. However, as we have already mentioned, in the next section we will 
give an example which shows that there may be a fuzzy simulation/bisimulation of a given type between 
two fuzzy automata, and there is not any crisp simulation/bisimulation of this type between them. 

It is worth noting that functions ((/)''T/ for all w e [fs, bs, fb, bb, fbb, bfb], can be characterized as follows: 

{a, b) e {(pfJio) « (Vx e X)(Vfl' e A) b^{a, a') ^ (6f o Q-^){b, a'), 

{a, b) e {(p^Jio) o (Vx e X)(Vfl' e A) b^{a' , a) < (p o 6f )(«', b), 

{cpffiQ) = (cpfJiQ) A [{cj}fj{g-')]-\ {(P^Yig) = (cp'rig) a l{(p'r{g-')]-\ 

i^^riQ) = ((ph'ig) A licp'rig-')]-', (cp'^'Yig) = i<P'Jig) a [{cpfj{g-')]-\ 

for all g e ??"(A, B),aeA and beB. 



6. Computational examples 

In this section we give examples which demonstrate the application of our algorithms and clarify rela- 
tionships between different types of simulations and bisimulations. 

The first example demonstrates the case when there are all types of simulations and bisimulations, but 
there is not any crisp bisimulation between two given automata. 

Example 6.1. Let = (A, 6^, 0^,1"^) and S = (B,6^,a^,T^) be automata over the Godel structure and an 
alphabet X - [x, y], with |A| = 3, \B\ - 2, and fuzzy transition relations and fuzzy sets of initial and terminal 
states which are represented by the following fuzzy matrices and vectors: 



o^ = [l 1 1], 



5^ = 



6' = 



1 

0.5 
0.4 

1 

0.6 



0.3 

1 

0.6 

0.6 
0.7 



0.4 
0.3 
0.7 



0.5 
0.3 
0.7 

0.6 
0.7 



0.6 
0.3 
0.7 

0.6 

1 



0.2 
0.4 
1 



T = 



(57) 
(58) 
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Using algorithms based on Theorem l5.3l we obtain that there are all types of simulations and bisimulations 
between fuzzy automata ^ and S, and the greatest simulations and bisimulations are given by matrices 



' 1 


0.7 




' 1 


0.7' 




' 1 


0.6' 




" 1 


0.7' 




" 1 


0.7' 




' 1 


0.6' 


1 


0.7 


, cp'"' = 


1 


0.7 




1 


0.6 


, <p'"' = 


1 


0.7 




1 


0.7 




1 


0.6 


0.6 


1 




0.7 


1 




0.6 


1 




0.7 


1 




0.6 


1 




0.7 


1 



On the other hand, using the version of the algorithms for crisp simulations and bisimulations, we obtain 
that there is not any crisp bisimulation between fuzzy automata ^ and S. 

The second example demonstrates the case when there are is a forward bisimulation, but there is not 
any other type of bisimulations between two given automata. 

Example 6.2. Let us change in l|57t and in l|58|l to 

a^ = [0 1 O], aS = [l 0.5]. 

Then the greatest forward bisimulation between fuzzy automata Jl and S is given by 



(p 



1 0.6 
1 0.6 
0.6 1 



but there is not any other type of bisimulations between J{ and S. 

Due to the duality between forward and backward bisimulations, it is possible to construct automata 
between which there is a backward bisimulation, but there is not any other type of bisimulations. 

We can also give an example which demonstrates the case when there are is a backward-forward bisim- 
ulation, but there is not any other type of bisimulations between two given automata. 

Example 6.3. Let us change in (|57t and in (|58t to 



[O l]. 



[0.7 l]. 



Then the greatest backward-forward bisimulation between fuzzy automata J{ and S is given by 



1 0.6 
1 0.6 
0.7 1 



but there is not any other type of bisimulations between J{ and S. 

Due to the duality between backward-forward and forward-backward bisimulations, it is possible to con- 
struct automata between which there is a backward bisimulation, but there is not any other type of bisim- 
ulations. 

Next, we give an example which demonstrates the case when there is not any type of simulations and 
bisimulations between two given automata. 

Example 6.4. Let us change in ll57t and in ll58ll to 

= O], ff^ = [0.5 l]. 

Then there is not any type of simulations and bisimulations between J\ and S. 

The following example considers the case (over the product structure) when the sequence of fuzzy rela- 
tions defined by (|49)l is infinite, and its intersection is the greatest forward bisimulation between the given 
fuzzy automata. 
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Example 6.5. Let J?I = (A, b^, a^, t^) and S = (B, 6^, a^, t^) be automata over the Goguen (product) structure 
and an alphabet X = {x}, with |A| = 3, \B\ - 2, and fuzzy transition relations and fuzzy sets of initial and 
terminal states which are represented by the following fuzzy matrices and vectors: 



[l 1 l] , 





1 


1 







1" 




1 


1 







1 










1 

2i 




1 



1 







1 





1 

2. 




1 



Computing the sequence l(pk]k€K for forward bisimulations by the formula ||49]| (w = fb) we obtain that 



1 ^1 

1 

^ 2k-l 

1 



, for each A: e N, 



= f\n = 



1 
1 
1 



According to Theorem l5.4[ cp is the greatest forward bisimulation between fuzzy automata J?I and S. 

The last example shows that the finiteness of the subalgebra (lm(!/'"')UlJ.v6x( Im(6^)Ulm(6f ) )) of X, which 
appears as an assumption in Theorem l5.3l is sufficient for the finiteness of the sequence defined by l|49] |, but 
it is not necessary. 



Example 6.6. Let us change a^, , and in the previous example to 

a^ = [l 1 O], a^ = [l O], 
Computing the fuzzy relations cp^, fc e N, using the formula (|49] |, we obtain that 



T = 



1 

1 
1 



and it is the greatest greatest forward bisimulation between fuzzy automata J{ and S. 

On the other hand, we have that the subalgebra (^(1/^'") U lm(6^) U lm(6f ) ) of X is not finite, since it 
contains ^, for every fc e N. 



7. Concluding remarks 

Simulations and bisimulations between fuzzy automata have been introduced and studied recently in fl^ . 
They have been defined as solutions of particular systems of fuzzy relation inequalities. Using the method 
for computing the greatest solutiuons of similar systems of fuzzy relation inequalities developed in [33], in 
the present paper we provide efficient algorithms for deciding whether there is a simulation or bisimulation 
of a given type between given fuzzy automata, and for computing the greatest simulation/bisimulation of 
this type, if it exists. 

Related systems of fuzzy relation inequalities will be discussed in the general context in our further 
research. 
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